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Solid  hydrogen  and  its  isotopomers  (HD  and  D2)  have 
long  been  recognized  as  unique  cryogenic  media  for  high- 
resolution  spectroscopic  studies  of  molecular  rovibrational 
dynamics  in  condensed  phases  [1,2].  In  the  lowest  energy 
state  of  the  solid  H2  crystal,  solid  parahydrogen  (PH2), 
each  H2  molecule  is  in  its  j=0  rotational  state  and  is  there¬ 
fore  a  spherically  symmetric  object  with  no  electrostatic 
multipole  moments.  Consequently,  the  pH2  crystal  is 
bound  together  only  by  weak  pH2~pH2  dispersion  interac¬ 
tions,  making  solid  pH2  a  very  “soft”  and  nearly  non¬ 
perturbing  environment  for  molecular  impurities.  This  has 
motivated  recent  high-resolution  infrared  (IR)  absorption 
studies  [24]  of  the  rovibrational  spectra  and  dynamics  of 
dopants  embedded  in  solid  pH2  matrices.  Analysis  of 
these  spectra  within  the  framework  of  crystal  field  theory 
provides  valuable  information  on  dopant-pH2  interactions 
and  on  the  microscopic  nature  of  dopant  trapping  sites  in 
the  pH2  matrix. 

Complementary  information  can  be  obtained  by  in¬ 
vestigating  the  impurity-induced  changes  in  the  IR  ab¬ 
sorption  spectrum  of  the  H2  matrix  itself.  For  instance, 
pH2  solids  containing  low  concentrations  of  j=l  orthohy¬ 
drogen  (0H2)  rotational  “impurities”  exhibit  an  IR  absorp¬ 
tion  feature  near  4153  cm'*  that  is  assigned  to  the  pure  vi¬ 
brational  Qi(0)  transition  (v=l  0,  j=0  ■<—  0)  of  pH2 
molecules  in  the  matrix  [5].  Although  this  transition  is 
strictly  IR  inactive  in  isolated  gas  phase  pH2  molecules 
and  in  pure  solid  pH2,  orientational  transitions  of  the  0H2 
dopant's  permanent  quadrupole  moment  provide  a  mecha¬ 
nism  for  inducing  nonzero  Qi(0)  transition  moments  in 
pH2  molecules  within  the  doped  solid  [6].  This  0H2- 
induced  absorption  feature  has  a  characteristic  asymmetric 
lineshape  arising  from  “vibron  hopping,”  or  delocalization 


of  the  v  =  1  vibrational  excitation  throughout  the  pH2  ma¬ 
trix.  Analysis  of  this  lineshape  provided  important  infor¬ 
mation  on  the  vibrational  dependence  of  the  H2— H2  poten¬ 
tial  before  the  era  of  ob  initio  and  molecular  beam  studies 
of  intermolecular  potentials  [6-8]. 

Spherical  dopants  such  as  alkali  metal  and  rare  gas 
atoms  can  .also  induce  IR  activity  in  solid  pH2  matrices  [9], 
although  via  a  qualitatively  different  mechanism  that 
originates  in  short-range  isotropic  overlap-induced  transi¬ 
tion  dipole  moments  [10]  arising  from  intermolecular  ex¬ 
change  interactions.  The  effectiveness  of  this  induction 
mechanism  stems  from  the  fact  that  a  chemical  impurity 
breaks  the  local  symmetry  of  the  pH2  lattice;  hence  this 
mechanism  is  truly  general,  and  applies  to  atomic  and 
molecular  impurities  alike.  Because  the  net  transition 
moment  for  these  induced  transitions  depends  sensitively 
on  the  microscopic  structure  of  the  dopant  trapping  site 
[10],  careful  study  of  these  spectral  features  can  provide 
new  insights  into  the  nature  of  impurity  solvation  in  the 
pH2  quantum  solid. 

In  this  Letter,  we  report  high-resolution  investigations 
of  the  IR  spectra  of  solid  pH2  matrices  doped  with  N2  im¬ 
purities.  Because  both  j=0  and  j=l  N2  rotational  levels  are 
populated  in  the  cryogenic  solid  pH2  environment,  these 
doped  matrices  turn  out  to  be  good  model  systems  for 
comparing  the  newly-identified  isotropic  induction 
mechanism  with  previously  studied  mechanisms  based  on 
transition  dipoles  induced  by  the  electrostatic  fields  of  im¬ 
purity  species.  This  is  because  the  j=0  N2  dopants  are 
spherical  objects  with  no  multipolar  electrostatic  field,  and 
thus  act  only  via  the  short-range  isotropic  mechanism  to 
induce  IR  activity  in  adjacent  H2  molecules,  while  the  j=l 
N2  dopants  (like  0H2  molecules)  induce  transition  dipoles 


in  both  nearby  and  distant  pH2  molecules  by  virtue  of  their 
permanent  quadrupole  moments.  Henc^two  distinct  in¬ 
duction  mechanisms  are  active  in  a  single  sample. 

We  prepare  millimeters-thick  doped  solid  pHi  sam¬ 
ples  by  rapid  vapor  deposition  [9]  of  flows  of  precooled  H2 
gas  and  room  temperature  dopant  gas  onto  a  BaF2  sub¬ 
strate  cooled  to  T  «  2  K  by  a  liquid  helium  bath  cryostat. 
The  precooled  H2  gas  emerges  from  a  variable  temperature 
ortho/para  H2  converter  [1 1]  typically  operated  at  T  «  15  K 
to  produce  nearly  pure  pH2  gas  flows  with  approximately 
100  parts  per  million  (ppm)  of  residual  0H2.  In  some  cases 
we  operate  the  converter  at  elevated  temperatures  to  pro¬ 
duce  oH2-enriched  samples. 

The  absorption  spectra  of  the  doped  pH2  samples  are 
recorded  along  the  substrate  normal  using  a  Fourier  trans¬ 
form  IR  spectrometer.  As-deposited  spectra  are  recorded 
immediately  after  deposition  at  a  substrate  temperature  of 
T  «  2.4  K.  We  then  anneal  the  samples  by  raising  the  sub¬ 
strate  temperature  to  T  «  4.8  K,  and  record  some  spectra  at 
this  temperature  as  well.  Sample  thicknesses  are  measured 
to  within  ±  10  %  (95  %  confidence  level)  using  the  ap¬ 
proach  described  in  Ref.  12.  The  impurity  concentrations 
cited  here  are  ratios  of  the  quantities  of  dopant  and  H2  en- 


FIG.  1.  Absorption  spectra  of  four  as-deposited  doped  pHi  sol¬ 
ids  at  T  «  2.4  K.  Trace  (a)  is  for  a  sample  containing  10  ppm  of 
CH4;  trace  (b)  is  for  a  sample  containing  10  ppm  of  CH4  and 
2500  ppm  of  0H2;  trace  (c)  is  for  a  sample  containing  10  ppm  of 
CH4  and  1000  ppm  of  N?;  trace  (d)  is  for  a  sample  containing 
10  ppm  of  CH4  and  1000  ppm  of  Ar.  Sample  thicknesses  are  2.7, 
.4.0, 3.1,  and  2.9  mm  for  traces  (a)  ftrough  (d)  respectively. 

Trace  (b)  has  been  scaled  by  a  multiplicative  factor  of  0.3  to 
compensate  for  the  higher  dopant  concentration  and  thickess  of 
this  sample.  Traces  (a)  throu^  (d)  are  presented  at  respective 
resolutions  of  0.0075,  0.008,  0.01,  and  0.008  cm'‘  and  have  been 
displaced  vertically  for  ease  of  presentation. 


tering  the  sample  chamber.  These  values  may  differ  from 
the  actual  impurity  concentrations  in  the  solid  (due  to 
varying  sticking  efficiencies  of  pH2  and  dopant  molecules, 
among  other  reasons)  but  are  estimated  to  be  within 
±  40  %  of  the  in  situ  concentrations. 

Figure  1  shows  the  ER  absorption  spectra  of  four  as- 
deposited  pH2  samples.  Trace  (a)  depicts  the  spectrum  of  a 
pH2  sample  containing  10  ppm  of  CH4,  which  is  included 
in  each  sample  as  a  structural  tracer  for  monitoring  the 
morphology  of  the  as-deposited  samples  [3,4].  No  detect¬ 
able  IR  activity  in  the  pHa  Qi(0)  region  is  induced  either 
by  the  CH4  dopant  (at  this  concentration)  or  by  the  residual 
0H2  molecules  present  in  the  sample. 

Trace  (b)  shows  the  absorption  spectrum  of  solid  pHa 
containing  2500  ppm  of  0H2,  and  displays  both  the  0H2- 
induced  pHa  Qi(0)  feature  between  4152  and  4153.2  cm  * 
and  the  narrow  Qi(l)  transition  of  the  dopants  themselves 
at  4146.6  cm*'.  The  Qi(0)  feature  exhibits  the  asymmetric 
vibron  hopping  lineshape  mentioned  previously.  This 
characteristic  lineshape  maps  out  the  subset  of  delocalized 
Qi(0)  vibrons  whose  spatial  wave  functions  achieye 
nonzero  overlap  with  the  quadrupolar  electrostatic  field  of 
the  oHa  dopant  [6],  and  therefore  indicates  the  presence  of 
an  induction  mechanism  based  on  quadrupole-induced 
transition  moments. 

Trace  (c)  depicts  the  spectrum  of  a  pHa  sample  con¬ 
taining  1000  ppm  of  natural  isotopic  abundance  NaJ  it 
shows  a  broad,  asymmetric  absorption  feature  from  about 
4150  to  4154  cm*'.  The  position  of  this  feature  coincides 
closely  with  that  of  the  pHa  Qi(0)  transition  observed  in 
mixed  oHa/pHa  samples,  and  we  assign  this  feature  to  Na- 
induced  Qi(0)  transitions  of  pHa  molecules  in  the  doped 
solid. 

The  lineshape  of  this  feature  shares  some  common 
aspects  with  the  oHa-induced  Qi(0)  feature  shown  in  trace 

(b) ,  including  a  relatively  steep  blue  edge  and  a  gently 
sloping  red  edge  that  meet  at  an  absorption  maximum  at 
4153.1  cm*'.  Despite  these  similarities,  we  can  rule  out 
any  possibility  that  the  feature  shown  in  trace  (c)  may  arise 
from  residual  oHa  impurities,  because  the  residual  oHa 
concentration  in  the  Na-dopsd  sample  is  too  low  to  pro¬ 
duce  the  Qi(l)  feature  seen  at  4146.6  cm*'  in  the  oHa- 
enriched  sample.  Furtheimore,  the  feature  shown  in  trace 

(c)  is  much  broader  at  the  baseline  than  is  the  oHa-induced 

Qi(0)  transition  shown  in  trace  (fr).  We  therefore  attribute 
the  Na-induced  feature  in  trace  (c)  to  the  superposition  of 
two  Na-induced  pHa  Qi(0)  transitions:  a  quad- 

rupole-generated  feature  induced  byj=l  Na  dopants,  which 
gives  the  central  portion  of  the  peak  its  characteristic 
shape,  and  a  broader,  unstructured  feature  induced  by 
spherical  j=0  Na  dopants. 

Support  for  this  interpretation  is  provided  by  trace  (d) 
in  Fig.  1,  which  depicts  the  IR  spectrum  of  a  pHa  sample 
containing  1000  ppm  of  Ar.  The  spherical  Ar  dopants  in¬ 
duce  a  broad,  structureless  pHa  Qi(0)  feature  visible  from 


4150  to  4153  cm*',  much  like  the  underlying  broad  feature 
in  trace  (c)  that  we  have  attributed  to  j=0  N2  molecules. 

As  we  explain  in  more  detail  below,  the  lineshapes  of 
spectral  features  attributed  to  the  isotropic  overlap  induc¬ 
tion  mechanism  differ  from  those  of  features  arising  from 
quadrupole  induction  because  different  subsets  of  the  pH2 
crystal's  Qi(0)  vibrons  are  rendered  IR  active  by  the  two 
mechanisms.  This  shows  that  careful  study  of  these 
dopant-induced  features  can  shed  new  light  on  the  funda¬ 
mental  physics  of  the  host  pH2  matrix. 

The  Qi(0)  features  we  have  attributed  to  the  symme¬ 
try  breaking  isotropic  induction  mechanism  originate  in 
(vibrational)  transition  dipole  moments  induced  by  short- 
range  intermolecular  exchange  interactions  [10].  The 
same  transition  moments  lead  to  collision-induced  Qi(0) 
H2  IR  absorption  in  pure  H2  [13,14]  and  in  gaseous  mix¬ 
tures  of  H2  with  other  atomic  and  molecular  species  [15], 
and  the  Qi(0)  pH2  features  induced  by  Ar  and  j=0  N2 
dopants  thus  represent  solid-phase  analogues  of  collision- 
induced  absorption.  The  “collisions”  which  generate 
overlap-induced  transition  moments  in  the  pH2  matrix  are 
simply  intimate,  short-range  interactions  between  nearest- 
neighbor  molecules  in  the  solid.  These  nearest-neighbor 
interactions  play  an  especially  prominent  role  in  solid  pH2 
matrices  because  the  individual  pH2  molecules  in  the  solid 
undergo  large-amplitude  zero  point  motions  about  their 
nominal  lattice  sites  [17];  this  in  turn  amplifies  the  effects 


FIG.  2.  Temperature  dependence  of  the  absorption  spectra  of 
the  3.1-mm-thick  N2/PH2  sample  shown  in  trace  (c)  of  Fig.  1. 
Trace  (a)  is  as-deposited  at  T  »  2.4  K;  trace  (b)  is  recorded  at  T  « 
4.8  K.  Both  spectra  are  presented  at  0.1  cm"'  resolution  and 
have  been  displaced  vertically  for  ease  of  presentation.  The  no¬ 
tations  +So(j)  and  +O0G)  represent  cooperative  transitions  in 
which  the  pHi  Qi(0)  vibrational  excitation  is  accompanied  by  a 
pure  rotational  transition  of  the  N2  dopant  from  the  (v  =  0,  j) 
state  to  the  (v  =  0,  j+2)  and  (v  =  0,  j-2)  states,  respectively. 


of  overlap-induced  transition  dipole  moments,  which  de¬ 
crease  rapidly  with  increasing  intermolecular  distance 
[16], 

Because  these  overlap-induced  Qi(0)  transition  mo¬ 
ments  depend  strongly  on  the  dopant“pH2  distance,  the 
transitions  rendered  IR  active  by  spherical  dopdnts  are 
those  in  which  the  final  state's  Qi(0)  vibron  has  substantial 
amplitude  on  one  or  more  of  the  twelve  pH2  molecules 
which  are  adjacent  to  the  dopant.  (In  addition,  the  final 
state's  vibron  wave  function  must  have  appropriate  sym¬ 
metry,  so  that  the  transition  moment  vectors  associated 
with  these  twelve  pH2  molecules  do  not  sum  to  zero.)  In 
contrast,  the  electrostatic  field  of  0H2  and  j=l  N2  dopants 
decays  rather  slowly  with  increasing  distance  (as  l/R'*), 
and  the  quadrupole  induction  mechanism  therefore  also 
activates  vibrons  with  substantial  amplitudes  on  non- 
nearest-neighbor  pH2  molecules.  Consequently,  the  spec¬ 
tral  features  induced  by  spherical  dopants  give  us  unique 
insight  into  the  perturbation  of  the  host  crystal's  vibron 
bands  by  chemical  impurities,  making  them  important 
tools  for  studying  solvation  in  the  pH2  matrix  environment. 

Our  interpretation  of  trace  (c)  in  Fig.  1  assumes  that 
N2  dopants  in  solid  pH2  retain  good  rotational  quantum 
numbers.  Evidence  that  this  is  so  is  provided  in  Fig.  2, 
which  shows  the  temperature  dependence  of  four  satellite 
IR  absorption  features  in  the  vicinity  of  the  N2-induced 
Qi(0)  feature.  The  strongest  of  these  features  is  a  multiplet 
that  decreases  in  intensity  upon  heating,  consisting  of  three 
peaks  at  4162.3, 4163.2,  and  4164.9  cm"',  each  0.9  cm"’ 
wide.  A  weak  structureless  feature  whose  intensity  is  in¬ 
dependent  of  temperature  appears  centered  at  4170.6  cm"' 
and  is  1 .3  cm"'  wide.  Two  very  weak  features,  each  about 
4.2  cm"'  wide,  appear  when  the  sample  is  heated  to 
T  »  4.8  K;  one  is  a  doublet  centered  at  4142.1  cm"'  and 
the  other  is  a  shoulder  centered  at  4178.5  cm"'  superim¬ 
posed  on  the  the  pH2  Qr  phonon  sideband.  The  spectral 
changes  produced  by  heating  the  doped  pH2  sample  are 
fully  reversible. 

As  we  show  next,  the  observed  temperature  depend¬ 
ence  of  these  satellite  features,  and  their  positions  relative 
to  the  sharp  maximum  of  the  Qi(0)  peak,  indicate  that  they 
arise  from  cooperative  transitions  in  which  the  pH2  Qi(0) 
transition  is  accompanied  by  a  pure  rotational  Aj  =  2  tran¬ 
sition  of  the  N2  dopant.  This  demonstrates  that  N2  dopants 
rotate  nearly  freely  in  the  pH2  matrix. 

The  homonuclear  N2  molecule  exists  in  ortho  and 
para  nuclear  spin  modifications  associated  with  even  and 
odd  rotational  quantum  number^>respectively  [18].  In  the  - 
absence  of  a  catalyst,  interconvefsion  between  orthonitro¬ 
gen  and  paranitrogen  is  very  slow;  hence  the  even  and  odd 
j  levels  ofN2  achieve  thermal  equilibrium  independently  in 
the  cryogenic  pH2  matrix.  Although  the  rotational  con¬ 
stant  of  N2  dopants  in  the  matrix  environment  differs 
slightly  from  that  of  isolated  gas  phase  N2  molecules  (see 
below),  this  difference  is  small  enough  that  we  can  esti- 


mate  the  rotational  level  populations  of  N2  dopants  rea¬ 
sonably  well  using  the  gas  phase  rotational  constant  Bgas  = 
1.99  cm-*  [19],  At  T  «  2.4  K  only  the  j=0  and  j=l  levels  of 
N2  are  appreciably  populated,  with  fewer  than  0.5  %  of  the 
N2  molecules  in  levels  with  j  ^  2.  Heating  the  sample  to 
4.8  K  promotes  approximately  15  %  of  the  j=0  molecules 
to  the  j=2  level,  but  excites  fewer-than  1  %  of  the  j=l 
molecules  to  the  j=3  level. 

Consequentlj^.the  triplet  feature  shown  in  Fig.  2, 
which  decreases  in'  intensity  by  roughly  20  %  upon  heat¬ 
ing,  can  be  assigned  to  cooperative  transitions  which  com¬ 
bine  the  Qi(0)  pH2  excitation  with  So(0)  Aj  =  2  transitions 
of  the  N2  dopant  from  its  j=0  level.  The  relative  intensities 
of  the  three  peaks  remain  constant  during  the  annealing 
cycle,  confirming  that  the  peaks  arise  from  a  common  ini¬ 
tial  state  that  is  depopulated  upon  heating.  The  triplet  na¬ 
ture  of  this  feature  suggests  that  interactions  between  the 
pH2  vibron  and  the  N2  dopant's  j=2  rotational  wave  func¬ 
tion  lift  the  orientational  degeneracy  of  the  j=2  state. 

Given  the  approximate  1 :2:2  intensity  ratio  of  the  three 
components  of  the  triplet,  we  assign  the  peaks  at  4162.3, 
4163.2,  and  4164.9  cm'*  to  transitions  to  the  mj  =  0,  ±  1, 
and  ±  2  levels  of  the  j=2  state,  respectively. 

The  weak  features  that  appear  only  in  the  spectrum  of 
the  heated  sample  can  be  assigned  to  cooperative  transi¬ 
tions  involving  Aj  =  ±  2  transitions  from  the  N2  j=2  level, 
as  indicated  by  the  Oo(2)  and  So(2)  labels  in  Fig.  2.  This 
N2  rotational  level  becomes  populated  (at  the  expense  of 
the  N2  j=0  level)  when  the  sample  is  heated  to  T  «  4.8  K. 
Conversely,  the  temperature-independent  feature  at  4170.6 
cm"'  arises  from  cooperative  transitions  involving  So(l)  Aj 
=  2  transitions  from  the  N2  j=l  level;  the  population  of  this 
Nt  level  is  unaffected  by  the  thermal  cycling  depicted  in 
Fig.  2. 

To  confirm  these  rotational  assignments,  we  note  that 
the  center  of  gravity  of  the  oH2-induced  Qi(0)  feature 
[trace  (b)  in  Fig.  1],  at  4152.8  cm"*,  marks  the  vibron-free 
origin  of  the  quadrupole-induced  Qi(0)  vibron  band  in 
solid  pH2  [6];  this  represents  the  energy  required  to  excite 
a  Qi(0)  transition  localized  on  a  single  pH2  molecule  in  the 
solid.  The  centers  of  gravity  of  the  N2-induced  satellite 
features  are  shifted  from  this  band  origin  by  amounts  cor¬ 
responding  to  the  energies  of  the  appropriate  N2  Aj  =  ±  2 
rotational  transitions,  which  are  simply  integer  multiples  of 
the  N2  rotational  constant  B  in  the  matrix  environment:  6 
B  for  the  So(0)  and  Oo(2)  transitions,  10  B  for  the  So(l) 
transition,  and  14  B  for  Ae  So(2)  transition.  Applying  this 
analysis  to  the  features  shown  in  Fig.  2  yields  B  =  1.79  ± 
0.04  cm"',  which  is  reduced  by  10  %  from  Bgas  1-99  cm'* 
by  crystal  field  effects.  The  good  agreement  among  the 
rotational  constants  derived  from  the  four  satellite  features 
indicates  that  the  rotational  assignments  in  Fig,  2  are  inter¬ 
nally  consistent. 

To  summarize,  we  have  presented  high-resolution  IR 
absorption  spectra  of  N2-  and  Ar-doped  solid  pH2.  The 


dopant-induced  features  observed  in  the  spectra  provide 
information  about  the  solvation  of  small  dopants  in  the 
pH2  matrix  and  about  the  symmetry  requirements  and 
length  scale  of  the  induction  mechanism  generating  the  IR 
activity.  In  particular,  the  lineshapes  of  the  dopant- 
induced  features  tell  us  which  vibrons  of  the  pH2  crystal 
are  rendered  active  by  a  particular  dopant,  and  how  the 
crystal's  vibron  bands  may  be  perturbed  by  the  dopant's 
presence. 

The  information  encoded  in  these  dopant-induced  IR 
absorption  spectra  can  be  extracted  with  the  help  of  quan¬ 
tum  Monte  Carlo  simulations  [10]  which  account  for  the 
highly  correlated,  large-amplitude  zero  point  motions  of 
the  pH2  solute  molecules.  Simulations  of  N2-doped  pH2 
should  enable  a  comprehensive  analysis  of  the  lineshapes 
and  multiplet  splittings  of  the  N2-induced  spectral  features, 
thereby  shedding  light  on  the  rotational  dynamics  and  sol¬ 
vation  of  chemical  impurities  in  a  highly  quantum  solid. 
Further  work  in  this  direction  is  in  progress. 
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